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COMPARATIVE  NUCLEOPHILIC  SUSCEPTIBILITY 
OF  ^.ELECTED  PHOSPHORO-  AND  PHOSPHONOTHIOLATES 
BY  SEMI-EMPIRICAL  MOLECULAR  ORBITAL  CALCULATIONS 


1 .  INTRODUCTION 
1 . 1  Problem  Definition. 

0-ethyl- S- (2 -diisopropylaminoethyl) -methylphosphonothiolate  (VX) 
(Figure  1)  is  a  potent  inhibitor  of  the  eni,yme  acetylcholinesterase  (AChE) . 
The  inhibition  mechanism  involves  phoshorylation  of  the  esteratic  site  in  the 
enzyme.  Because  of  the  lethal  toxicity  of  VX,  severe  restrictions  have  been 
placed  on  its  use  in  laboratory  experimentation. 


Figure  1.  VX 


In  an  effort  to  circumvent  many  difficulties  involving  the  use  of 
chemical  agents  in  the  laboratory,  the  U.S.  Army  Chemical  Research, 

Development  and  Engineering  Center  (CRDEC)  has  been  the  lead  organization  for 
the  development  and  use  of  chemical  agent  simulants  (conversation  with  William 
White,  Research  Directorate,  CRDEC,  May  1990).  Simulants  are  chemical 
compounds  that  possess  similar  properties  and/or  activities  as  an  agent  but 
are  non- toxic.  In  this  way,  simulants  will  behave  similarly  to  agents  in  a 
given  experiment  but  require  none  of  the  precautions  associated  with  the  use 
of  toxic  chemicals.  To  promote  the  use  of  simulants  and  provide  a  centralized 
location  to  find  properties  regarding  simulants,  CRDEC  created,  as  well  as 
maintains,  the  Chemical  Agent  Simulant  Data  Center.’ 

One  area  with  little  success  in  developing  simulants  has  been  in 
simulating  reactivity,  either  oxidative  or  hydrolytic.  This  study  will 
present  a  new  approach  to  examining  this  problem. 
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1 . 2  Potential  Simulants. 


Figures  2-8  display  the  potential  simulants  that  have  been 
examined  in  this  study.  The  compounds  examined  with  the  exception  of  0- 
ethyl-S-ethyl  methylphosphonothiolate  (EA-5533)  (Figure  5)  are  common 
pesticides.*  Research  Directorate,  CRDEC,  (conversation  with  Leonard 
Szafraniec,  Research  Directorate,  CRDEC,  May  1990),  synthesized  jlA-5533. 

1 . 3  Nucleophilic  Displacement  In  VX . 

VX  detoxification  can  cause  a  number  of  different  mechanisms. 
This  study  will  focus  only  on  the  nucleophilic  displacement  of  the  alkylthio 
group  by  a  nucleophilic  attack  at  the  phosphorus.  The  products  obtained  by 
this  attack  are  a  function  of  the  strength  of  the  nucleophile,  pH,  and 
solvent.  The  VX  reaction  with  DS2  (sodium  hydroxide,  methoxyethanol , 
diethylenetriamine)  involves  the  substitution  of  the  2- (bis - isopropylamino) 
ethylthiol  group  by  the  2-methoxyethoxide  (OMC)  (equation  I')  .  In 
water/alcohol  solutions,  competitive  substitution  by  the  alcohol  can  occur 
(equation  2) 
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1 .4  Theoretical  Chemistry  Applications. 

Computational  chemistry  is  extensively  used  in  the  examination  of 
chemical  compound  properties.  Although  a  number  of  computational  methods  have 
been  used  successfully  to  predict  physical  and  biological  properties, 
prediction  of  chemical  reactivity  has  been  less  successful Theoretical 
methods  allow  for  the  facile  computation  of  numerous  electronic  properties  of 
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Figure  2.  Simulant  1,  0-Methyl -S-Methyl 
•  Aminophosphorothiolate 


Simulant  2,  0-Methyl -S-Methyl-  Figure  4.  Simulant  3,  0-Ethyl -0-ChlorO' 
N-Acetyl  Amidophosphorothiolate  Bromophenyl-S-Propyl 

Phosphorothiolate 


Figure  5.  Simulant  4,  0-Ethyl -S-Ethyl 
Methyl phosphonothi ol ate 


Figure  6.  Simulant  5,  Bis(0-Ethyl )- 
S-Ethyl  Phosphorothiolate 


Figure  7.  Simulant  6,  Bis  (0-Methyl )-  r  c,-™  ^  ^  ^ 

Chlorophenyl  Phosphorothiolate  Simulant  7.  Tnethyi  Phosphate 


molecules  that  can  be  used  to  infer  reactivity  and  used  to  further  compute 
reaction  related  properties.  In  effect,  the  ground  state  electronic  and 
steric  attributes  are  used  to  qualitatively  predict  how  well  a  reaction  can 
occur . 


One  difficulty  in  the  application  of  molecular  orbital 
calculations  is  the  sheer  amount  of  electronic  data  that  can  be  generated. 

Two  approaches  to  apply  molecular  orbital  calculations  have  been  examined  in 
this  study.  The  first  approach  uses  hierarchical  clustering  techniques  to 
find  the  compound  that  best  matches  those  electronic  properties  expected  to 
influence  the  nucleophilic  displacement.  In  the  second  approach,  molecular 
graphics  have  been  applied  to  condense  the  plethora  of  numerical  data  into  an 
easier  to  comprehend  description  of  the  electrostatic  potential.  The 
electrostatic  potential  has  been  used  successfully  in  this  laboratory  and  by 
others  to  qualitatively  assess  the  reactivity  of  a  substrate 

2  .  COMPUTATIONAL  PROCEDURE  . 

All  geometry  optimizations  and  calculations  of  electronic 
properties  have  been  performed  using  the  MNDO  algorithm  as  contained  in  MOPAC 
v5.0  on  a  Stardent  Titan  super  workstation.**  Graphical  display  of  the 
electrostatic  potentials  and  molecular  structures  were  done  using  the 
Molecular  Modeling  Analysis  and  Display  System  (MMADS)  program.’*  The  pattern 
recognition  program,  ARTHUR,  was  used  to  perform  all  hierarchical 
clustering . ” 


3.  RESULTS  AND  DISCUSSION 

Two  methods  were  employed  to  determine  if  theoretical  and 
computational  methods  could  be  useful  in  selecting  simulants.  If  the 
important  properties  that  determine  the  reactivity  of  VX  with  respect  to 
nucleophilic  displacement  could  be  adequately  matched  by  a  potential  simulant, 
then  the  reactivities  would  likely  be  similar.  However,  the  important 
properties  in  hydrolysis  were  not  intuitively  clear.  Seven  theoretically 
determined  electronic  properties  were  calculated  and  used  to  attempt  to  model 
reactivity  (described  in  section  3.1).  Using  pattern  recognition  methods,  it 
is  possible  to  determine  which  of  the  potent  simulants  best  match  all  or  some 
of  the  electronic  properties  of  the  agent.  A  second  approach  using  the 
molecular  electrostatic  potential  (MEP)  was  also  employed.  In  particular,  the 
difference  between  the  MEP  of  the  simulant  and  VX  was  examined.  The 
reactivities  will  more  likely  be  similar  as  the  MEP  map  approach  exact 
overlap . 


♦Famini,  George  R.,  "A  Computational  Comparison  of  the  Electronic 
Properties  of  G  Agents  with  Selected  Reaction  Simulants,"  In  Proceedings 
q£_  the^gSJ  Scientific  Conference  on  Chemical  Defense  Research.  CRDEC- 
SP-024,  U.S.  Army  Chemical  Research,  Development  and  Engineering  Center, 
Aberdeen  Proving  Ground,  MD,  September  1990,  unpublished  data. 
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3 . 1  Hierarchical  Clustering. 


The  electronic  properties,  calculated  by  the  MNDO  algorithm, 
expected  to  contribute  to  the  electrophilic  properties  in  VX  are  as  follows: 

•  heat  of  formation, 

•  the  energy  of  the  highest  occupied  molecular  orbital, 

•  the  energy  of  the  lowest  unoccupied  molecular  orbital , 

•  the  electronic  charge  of  the  phosphorus , 

•  the  electronic  charge  of  the  sulfur  atom, 

•  the  electronic  charge  of  the  P  -  0  oxgyen,  and 

•  the  dipole  moment. 

These  electronic  properties,  the  values  for  VX,  and  the  seven 
simulants  are  listed  in  Table  1. 

Three  dendrograms  were  produced  by  the  pattern  recognition  package 
ARTHUR.  The  first  dendrogram  (Figure  9)  compares  all  seven  electronic 
properties  for  VX  and  the  simulants.  The  dendrogram  shows  that  Simulant  4,  0- 
ethyl-S-ethyl  methylphosphonothiolate  is  the  best  fit  for  VX.  Simulants  3,  6, 
5,  and  1  are  the  second  best  choices;  Simulant  2  is  the  third  best  choice.  A 
second  dendrogram  (Figure  10)  compares  six  (the  heat  of  formation  was  not 
considered)  of  the  Modified  Neglect  of  Diatomic  Differential  Overlap  (MNDO) 
electronic  properties  for  VX  and  the  simulants. 

Simulants  1,  3,  4,  5,  and  6  were  the  first  choice  for  the  best  fit 
and  Simulant  2  the  second  choice.  A  third  dendrogram  compares  only  the  energy 
of  the  orbit  (ELUMO) ,  because  the  Frontier  Molecular  Orbital  Theory  predicts 
that  this  property  is  the  most  significant  in  Sn2  chemistry.  For  the  ELUMO 
dendrogram  (Figure  11),  the  ranking  is  Simulant  4  followed  by  Simulant  1  for 
the  best  fit  for  VX. 

3 . 2  Electrostatic  Potentials. 

Electrostatic  potentials  have  been  shown  to  be  a  good  indicator  of 
relative  reactivity  due  to  attraction  and  repulsion  of  point  charges.  The 
Molecular  electrostatic  Potential  (MEPs)  were  calculated  using  a  modified 
coulombic  expression  as  shown  in  equation  3. 

q. 

Ra-  r 


Vapp  =  ^ 
A 


V^  is  the  approximate  electrostatic  potential 
q,  is  the  formal  charge  on  atom  A 
Ry^  is  the  position  of  the  atom 
r  is  the  electron  at  point  r 


Electrostatic  potentials  for  VX  and  the  simulants  are  shown  in 
Figures  12-19.  The  planes  displayed  are  the  phosphorus,  the  double  bonded 
oxygen,  and  the  leaving  group.  The  significant  region  of  the  plot  is  the  area 
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Table  1.  MNDO  Electronic  Properties 


V£> 

o 

o 

to 

CO 

CO 

o 

r' 

o 

VO 

to 

o 

00 

rH 

s 

V£) 

(N 

H 

VO 

CO 

o 

OV 

Q 

• 

n 

• 

fO 

• 

VO 

• 

rH 

• 

rH 

• 

CO 

• 

rH 

CM 

o 

to 

00 

o 

r* 

00 

Vi 

00 

to 

VO 

o 

OV 

CO 

CM 

OV 

CT 

00 

fO 

rH 

00 

OV 

OV 

rH 

o 

rH 

o 

o 

o 

o 

1 

1 

1 

1 

r 

0) 

XJ 

00 

o 

to 

o 

rH 

o 

Its 

O 

o> 

rH 

CM 

<VJ 

CM 

OV 

VO 

o 

rH 

cr 

00 

o 

to 

CO 

CO 

o 

CO 

o 

Ifl 

• 

VO 

• 

to 

• 

VO 

• 

VO 

• 

VO 

• 

VO 

• 

VO 

• 

0)  -H 

■p  x: 

1 

1 

1 

1 

1 

• 

1 

' 

Its  -P 

fH  O 

rH 

0  P 

rH  (0 

•HO  0) 

lO  4-> 

x:  x;  4J 

n 

o 

to 

CO 

o 

CM 

-H 

-pa  Its 

cu 

rH 

r- 

00 

to 

CM 

OV 

\D 

•H  JQ 

0  U  rH 

tr 

ov 

OV 

O' 

rH 

00 

a\ 

U 

0)  P  O  0 

00 

ov 

o 

o 

00 

o 

o 

CO 

U  O 

-p  0  x:  fl)  -H 

• 

• 

• 

• 

« 

• 

• 

• 

O 

Its  x:  a-p  0)  x: 

rH 

rH 

rH 

rH 

rH 

rH  a  Its  P  -P 

^  (0 

O  W  iH  rH  10  o 

(0  rH 

•H  0  >1  0  rH  P 

■H  3 

x;  x:  a-p  o  o 

3  0 

p  a  0  x:  -p  x: 

fH 

CM 

OV 

CM 

VO 

rH 

0  0) 

0  0  p  P  x:  a 

0 

VO 

<N 

ov 

C' 

to 

to 

(U  rH 

p  TS  a  o  p  u 

e 

OV 

VO 

to 

ov 

rH 

to 

n 

rH  0 

o-p  1  coo 

0 

n 

n 

to 

o 

CO 

CO 

r> 

CO 

0  E 

x:  E  to  0  p  x: 

rH 

r* 

cn 

H 

CO 

VO 

o 

CM 

00 

E 

a  Its  1  x:  0  a 

u 

« 

• 

• 

• 

• 

• 

• 

• 

•0 

to  rH  ax: 

iri 

iH 

<N 

«N 

iH 

CM 

CM 

o 

•a  0) 

0  fh  >,  w  a-H 

1 

1 

1 

1 

1 

1 

1 

1 

dJ-H  E 

x:  ><  c  0  to  >1 

•H  a  o 

ap  0)  x:  o  c 

a  3  -p 

o  V  js  ax3  « 

3  O  (0 

c  w  arH  ax; 

0  U  EE 

•p  10  o  >1  a 

VO 

CO 

CO 

00 

t> 

to 

CO 

U  0  U)  0  0 

El  E  X:  rH  0 

r» 

CO 

(JV 

00 

VO 

VO 

CO 

o  C  3  -P  4J 

10  B  O  P  >,  P 

0 

o 

CO 

rH 

o 

iH 

rH 

OV 

00 

3  P  lO  IQ 

1  p  0)  x:  o  0) 

e 

o 

rH 

TT 

o 

o 

o 

rH 

o 

tp  iH  X3  E  P  rP  P 

0 

n 

tT 

H 

to 

CO 

rH 

0)  £  C  P 

>1  >1 0  v  x:  <0 

» 

« 

9 

9 

9 

• 

9 

• 

w  (0  a  v  3 

x; x:  p  .p  1  ox: 

<n 

o 

O 

C\ 

O 

o 

OV 

rH 

£  «  to  pvtM 

p  p  o  >v(o  1  a 

1 

rH 

rH 

1 

H 

rH 

1 

rH 

C  cn  >  o  >i<H 

0)  0)  rH  x;  1  —  tn 

1 

1 

1 

1 

1 

0  -H  0  XS  X  3 

E  E  X:  P  --rH  0 

•H  x:  iH  a  o  to 

I  1  U  CS  rH  >1X3 

4J 

to  (0  1  1  >ix:  a 

Its  0)  0)  V  a>  0)  -P 

1  1  o  to  x:  p 

E  x:  x:  x:  x:  x:  c 

rH  iP  1  1  p  0)  rP 

p  -P  -P  -P  -P  -P  0) 

>1  >,rH  iP  w  E  >1 

O  E 

x3  x:  >1  >1 1  1  x: 

o 

VC 

CM 

OV 

rH 

c^ 

o 

<M  <M  <M  C  C  C  0 

p  p  x3  x:  o  o  p 

to 

r* 

VO 

«-( 

00 

to 

00 

rH 

o  o  o  o  o  s 

«  «  p  P  ww  o 

O) 

CN 

rH 

00 

VO 

00 

VO 

IH 

£  X  u  u  to  to  -P 

X 

r* 

M 

00 

VO 

o 

r-l 

oo 

0  ^  ^  0)  0)  Q)  0) 
O'  O'  O'  O' 

1  1  1  1  -H-P  p 

CM 

O 

CM 

r' 

cn 

o 

o  o  o  o  n  oQ  H 

• 

• 

• 

• 

• 

• 

• 

• 

^  ^  V4  ^4  0 

CO 

o 

CJV 

r» 

n 

to 

OV 

Its  0)  v  Its  Its  Its  a 

00 

ov 

OJ 

CM 

o 

to 

rH 

o 

0)  c  c  x:  x:  x:  -p 

P  CM  n  ^  to  VO  p- 

1 

1 

rH 

rH 

•H 

rH 

rH 

CM 

B  H  u  u  o  u  o 

1 

1 

' 

1 

1 

1 

II  II  11  II  II  II  II 

p  p  p  p  p  p  p 

C  C  C  C  C  C  C 

10  10  ns  lO  lO  lO  10 

0  0 

rH  rH  rH  rH  fH  ^  rH 

X( 

E  E 

3  3  3  3  3  3  3 

c 

O  3 

E  E  B  E  E  E  E 

3 

x;  <-i  a  o  w  X 

B  H  H  cr  D*  cro 

•H  •H  ‘H  •H  -H  “H 

0 

a 

e 

0 

o 

VX 

rH 

fV) 

CO 

to 

VO 

CO  to  CO  to  CO  (0  to 

15 


*  « 

m 

m 

m 

* 

* 


m 

♦ 


8;: 

-  § 


0"  ♦  ♦ 
.  I  • 


«  *  « 
« 

* 

« 

* 

* 

« 

*  *  ♦  4t 


«  «  «  «  * 
i» 

* 

m 

9 

m 

9 

♦ 


9 

9 

9 

9 

9 

9  9  9  9 


o  ^ 
o 


^  o  o  o  o 

o  o  o  o 


8  0  0  0  - 

o  5  o 


~  - 

- 

v« 

^  — 

ID 

CO 

<o 

CM 

K 

v> 

to 

<0 

CM 

K 

c 

s 

«i« 

C 

c 

c 

«« 

c 

« 

3 

«« 

c 

« 

3 

c 

« 

3 

c 

C 

« 

C 

« 

3 

E 

• 

3 

E 

«S 

3 

E 

« 

3 

E 

m 

3 

E 

«B 

3 

m 

3 

F 

s 

E 

E 

X  E 

E 

1 

(7$ 

05 

o5 

S  05 

05 

C5 

<A 

05 

05 

16 


Figure  10.  Dendrogram  for  Six  Electronic  Propertle 
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Figure  16.  Electrostatic  Map/Simulant  4 


Figure  17 


Electrostatic  Map/Simulant  5 


Figure  18.  Electrostatic  Map/Siroulant  6 


Figure  19.  Electrostatic  Map/Simulant  7 


180°  to  the  leaving  group.  This  significant  region  is  shown  as  solid  contour 
lines  and  represents  a  positive  area  attractive  to  a  negative  charge. 

Simulants  that  best  fit  the  reactivity  of  VX  will  have  a  positive  region 
similar  in  shape  and  size. 

Computer  generated  electrostatic  potential  difference  maps  have 
been  made  between  VX  and  each  of  the  simulants,  thus  simplifying  the 
comparison.  The  simulants  that  have  the  best  electrostatic  overlaps  will 
result  in  difference  maps  with  the  smallest  positive  lobe  opposition  to  the 
leaving  group.  The  electrostatic  overlap  maps  are  shown  in  Figures  20-26. 

Examination  of  the  diffence  map  (Figure  23)  shows  that  Simulant  4, 
0-ethyl-S-ethyl  methylphosphonothiolate ,  would  closely  resemble  VX'in 
nucleophilic  reactivity.  However,  Simulant  1,  0-methyl-S-methyl 
aminophosphorothiolate ,  and  Simulant  5,  0- ethyl -0- ethyl  S-ethyl 
phosphorothiolate  (Figures  20  and  24,  respectively),  may  provide  reasonable 
reactivity. 


Simulant  4  (EA-5533)  and  Simulant  5  were  compared  (Figure  27)  to 
determine  the  difference  in  phosphorus  reactivity  between  a  phosphonothiolate 
and  a  similar  phosphorothiolate.  The  large  positive  lobe  shown  indicates  a 
large  difference  in  reactivity  between  the  two  phosphorus. 


Figure  27.  Difference  Map/Simulant  4  and  Simulant  5 


4.  SUMMARY 


Hieratchical  clustering  and  electrostatic  potential  techniques 
were  used  to  predict  the  similarity  in  reactivity  between  several  potential 
phosphorothiolate  simulants  and  VX.  The  technique  considered  the  Sn2 
displacement  of  a  leaving  group  attached  to  phosphorus.  Two  compounds,  0- 
methyl-S -methyl  aminophosphorothiolate  and  0-ethyl-S-ethyl 
methylphosphonothiolate  were  identified  as  suitable  simulants.  Table  2 
summarizes  the  ranking  of  simulants  for  HEP  and  hierarchical  clustering. 


20 


Figure  20.  Difference  Map/VX 
and  Simulant  1 


Figure  21.  Difference  Map/VX 
and  Simulant  2 


Figure  22.  Difference  Map/VX 
and  Simulant  3 
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Figure  23.  Difference  Map/VX 
and  Simulant  4 


Difference  Map/VX  Figure  25.  Difference  Map/VX 

and  Simulant  5  and  Simulant  6 
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Figure  26.  Difference  Map/VX 
and  Simulant  7 


Table  2.  Simulant  Ranking 


Heirarchical  Clustering 


Rank 


7  Electrostatic  6  Electrostatic  1  Electrostatic  MEP 

Properties  Properties  Property 


1  simulant  U 


simulants  1,3,  simulant  4  simulant  4 

4.5.6 


2  simula''.ts  1,3,  simulant  2 

5.6 


simulants 
1.2, 3, 5. 6 


simulant  1 
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